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Abstract 

We give an FPRAS for Holant problems with parity constraints and not-all-equal con- 
straints, a generalisation of the problem of counting sink-free-orientations. The approach 
combines a sampler for near-assignments of "windable" functions - using the cycle- unwinding 
canonical paths technique of Jerrum and Sinclair - with a bound on the weight of near- 
assignments. The proof generalises to a larger class of Holant problems; we characterise this 
class and show that it cannot be extended by expressibility reductions. 

We then ask whether windability is equivalent to expressibility by matchings circuits (an 
analogue of matchgates), and give a positive answer for functions of arity three. 

1 Introduction 

In this paper we will show that the following problem has an FPRAS (a type of approxima- 
tion algorithm - see Section . 

Name #ParityNAE 

Instance A multigraph G in which each vertex is labelled Even, Odd, or NAE 
Output The number of subsets F C E(G) such that: 

• each Even vertex has an even number of incident edges in F 

• each Odd vertex has an odd number of incident edges in F 

• each NAE vertex has at least one incident edge in F and at least one incident 
edge in E(G) \ F 

Theorem 1. There is an FPRAS for #ParityNAE. 



1.1 Relationships with other counting problems 

A sink-free orientation of a graph is a choice of orientation of each edge such that no vertex 
has out-degree zero. The problem #SFO is: given a graph, count the number of sink-free 
orientations. We can also allow "skew" edges, where the ends of the edge must both be 
oriented outwards or both oriented inwards. 

Bubley and Dyer studied #SFO and gave an FPRAS @]. They showed as a corollary 
that there is an FPRAS for counting solutions to a formula in conjunctive normal form in 
which every variable appears at most twice, which they showed is a #P-hard problem. The 
first part of their argument was a standard reduction to sampling - finding a fully polyno- 
mial almost uniform sampler (FPAUS) for sink-free orientations. Then, they constructed 
a Markov chain that converges to the uniform distribution on sink-free orientations, and 
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Figure 1: Reduction from #SFO to ^ParityNAE. The edge with two arrows is a skew edge. A 
sink-free orientation is illustrated with the corresponding set F draw in thick grey. 

bounded its mixing time using a two-stage path coupling argument. Colin, Pemantle and 
Propp later gave an exact sampler with 0(|V| • |F|) mean running time using a kind of 
rejection sampling [5]. 

A simple reduction from #SFO to #ParityNAE is illustrated in Figure [TJ showing that 
^ParityNAE generalises the problem of counting sink-free orientations in a graph (while also 
allowing parity constraints). Given an instance G of #SFO, label all the vertices NAE, 
subdivide each non-skew edge uv, label the new vertices (which we will refer to as "m u „" ) 
Odd, then attach a degree-one Odd vertex to each NAE vertex. This gives an instance G' 
of #ParityNAE. For all orientations O of G define a set Fo C E(G') by taking all edges 
attached to degree-one Odd vertices and all "heads" : for non-skew edge uv take um uv g F 
if and only if uv is directed towards u, and for skew edges uv take uv € F if and only if uv 
is oriented outwards. Each degree- two Odd vertex in G 1 has exactly one incident edge in 
Fo, and each NAE vertex in Fo has at least one incident edge in F, and if O is sink- free 
then each NAE vertex in Fo has at least one incident edge not in F. Furthermore, any 
F C F(G") satisfying these conditions is Fo for some sink-free orientation O. The function 
O i y Fo therefore gives a bijection from sink-free orientations of G to the set of subsets of 
F(G') that get counted by #ParityNAE. 

^ParityNAE, at least when restricted to bounded-degree graphs, is a type of Boolean 
Holant problem. Holant problems are a quite general type of graphical counting prob- 
lem. The constraints Odd, Even and NAE in ^ParityNAE are generalised to functions 
F: {0, l} k —> C, called (Boolean) signatures in this context. Relations R C {0, l} fe can 
also be used by considering the function R: {0, l} k — > {0, 1} that takes the value 1 exactly 
on elements of R. In this discussion we will take the codomain to be the set of complex 
numbers, but afterwards we will restrict to non-negative rational-valued signatures, which 
we call weight-functions. 

A Holant instance is a graph G equipped with a function F v : {0, 1} J " — > C for each 
vertex v, where J v is the set of edges incident to v. In fact we will want to allow self-loops, 
which calls for a slightly more complicated definition - see Section [3J We are interested in 
the total weight 

n f »m-o- 

xe{o,i} B < G > vev(G) 

For example, if all F v take the value 1 on vectors of Hamming weight 1, and take the value 
otherwise, then the total weight is just the number of perfect matchings of G, because a 
vector x G {0, 1} B ( G ) is the characteristic vector of a perfect matching of G if and only if 

Uvev(G) F vi x \jv) = !■ 

Given a finite set F of signatures, Holant(J r ) is the problem of evaluating the total weight 
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where G is given as input, and where we require that each F v is a copy of some F £ F: for 
some enumeration vi, . . . ,v k of J v we have F v (x) = F(x.(vi), . . . , x(wfe)) for all x £ {0, l} Jl '. 
For example, let F be the function defined by F(0,0, 1) = F(0, 1,0) = F(1,0,0) = 1 and 
F(i,j, k) = elsewhere. Then Holant({F}) is the problem of counting perfect matchings in 
degree-three graphs. 

For all positive integers k define Even^, Oddfc, NAEj. : {0, l} fc — > {0,1} by setting 
Everife(a;i, . . . , Xk) to be 1 if and only if x\ + - ■ -+Xk is even, setting Odd/ £ (xi, . . . , x k ) to be 1 
if and only if x±+- • -+Xk is odd, and setting NAE/ £ (xi, . . . , Xk) to be 1 if and only if 1 < x± + 
■ ■ -+Xk < k—1. The restriction of #ParityN AE to graphs of maximum degree at most d is then 
equivalent to Holant (.Fd) where Fd = {Eveni, Oddi, NAEi, . . . , Even^, Odd^, NAE<j}. 
By Theorem [TJ this problem has an FPRAS for each d. 

#ParityNAE is also a counting constraint satisfaction problem, at least when restricted to 
bounded-degree graphs. Roughly speaking, a instance of a constraint satisfaction problem 
is a list of constraints like "x V y, y A z, y = z", and we are interested in the number of 
configurations satisfying the constraints. In particular [5], for any finite set of signatures F, 
an instance of #CSP(F) is a set of variables V and a list of formal function applications 

Fi(v M , . . . ,«i, fcl ), . . . ,F s (v Sj i, . . . ,v s ,k s ) 

where Fi : {0, l} ki — > C is a function in F for each 1 < i < s, and Vij £ V is a variable for 
each 1 < i < s and 1< j < fc,. The value of this instance is 

m 

xG{0,l} v 1=1 

For information about the complexity of approximately evaluating #CSPs, see [S]. 

A ^=CSP(J-") instance can be drawn as a "dual constraint hypergraph" , which has vertices 
for each of the constraints c\,...,c s , and a hyperedge {c, | v = Xi.j} for each variable v 
(ignoring multiplicities for now). The dual constraint hypergraph is a graph if and only if 
every variable appears exactly twice. In this way, Holant(J r ) is the restriction of #CSP(J-") 
to read-twice instances. Note that the variables of the #CSP are the edges of the Holant 
instance; sometimes a #CSP is described in the opposite way, as the primal constraint 
hypergraph, with variables as vertices and s edges or hyperedges. 

We now recall the relationships between #CSPs and Holant problems in the Hadamard 
basis as discussed in [TH] . Note that while these equivalences are usually stated in the context 
of exact evaluation, the reductions just involve preprocessing the input, and so also apply 
in the context of approximate counting. Firstly, equality constraints can be used to break 
the read-twice restriction. Letting = 3 denote the function {0, l} 3 — > C taking the value 
1 on (0,0,0) and (1, 1, 1), and taking the value elsewhere, if =3 is in F then Holant(J r ) 
is equivalent to #CSP(J r ) 7, Proposition 5]. Secondly, let F: {0, l} fc — > C denote the 
Hadamard transform, defined by 

F(x u ...,x k ) = 2- k / 2 F(y 1 ,...,y k )(-ir^+-+ x *y». 
ye{o,i} fc 

Holant(J r ) is always equivalent to Holant({F | F £ F}); see Proposition 1] or [22] . Also, 

F = F for any F. So if F contains =5, then Holant (.F) is equivalent to #CSP({-F | F £ F}). 
But =3 is just Everi3 multiplied by a factor of \pl (which can be easily accounted for). 

Taking F to be the set Fd defined above, with d > 3, we find that the restriction of 
#ParityNAE to instances of degree at most d is equivalent to #CSP({F | F £ Fd})- By 
Theorem [T] this problem has an FPRAS for each d. In this sense, #ParityNAE generalises 
#SFO to a #CSP. Note that 6dd^(0) = 1/V5 and Odd^(O) = -I/a/2- So we get a class 
of FPRASes for #CSPs using functions with mixed signs. 

1.2 Techniques 

Like Bublcy and Dyer we will use Markov chains, but to bound the mixing time we will 
instead apply the canonical paths technique. More precisely, we will use a multicommodity 
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flow with cycle-unwinding as used by Jerrum and Sinclair 19 . They proved the following 
relevant result: for any polynomial p we can sample efficiently from the uniform distribution 
of perfect matchings, in graphs G satisfying 



number of matchings of order 4|F(G)| — 1 
number of matchings of order i|V(G)| 



<P(\V(G)\). (1) 



Recall that a matching of a graph is a set of edges not sharing any vertices, and a 
matching is perfect if it has order |U(G)|/2. A perfect matching is a satisfying assignment 
to a certain system of constraints: each edge is either IN or OUT, and every variable 
enforces a perfect matchings constraint, that exactly one of its incident edges is IN. From 
this perspective a natural question is: what weight-functions can we use instead of perfect 
matchings constraints? We show that Jerrum and Sinclair's result generalises in a certain 
sense to windable functions, defined as follows. 

Definition 2. For any finite set J and any configuration x £ {0, 1} J define to be the 
set of partitions of {i \ Xi = 1} into pairs and singletons. A function F : {0, 1} J Q-° 
is windable if there exist values B(x,y, M) > for all x, y £ {0, 1} J and all M £ .M^ey 
satisfying: 

1. F(x)F(y) = Emgm;^ y= M ) for all x, y e {0, 1} J , and 

2. B(x,y,M) = B(x©S,y9S,M) for allx 7 y £ {0,1} J and all S £ M £ M' x(By . 

Here x © S denotes the vector obtained by changing Xi to 1 — Xi for the one or two 
elements i in S. 

The next question is: what kinds of constraints guarantee a bound like ([I])? We give one 
answer: strictly terraced functions. 

Definition 3. A function F : {0, 1} J — > Q-° is strictly terraced if 

F(x) = ==> F(x © e. t ) = F(x © e^) for all x £ {0, 1} J and all i,j £ J. 

Here x © ej denotes the vector obtained by changing xi to 1 — xi . 

We will discuss these definitions more throughout the paper. Using properties of these 
classes, we will establish Theorem [TJ A feature of the techniques is that they cannot be 
extended by expressibility reductions, where we just substitute a constraint by a "circuit" , 
a gadget gluing together other constraints. The following theorem makes this precise. 

Theorem 4. Let T be the class of strictly terraced windable functions. Then 

• T is closed under taking weight-functions of connected circuits 

• T contains Even^, Odd^, and NAEj for all k > 1 

• for all finite subsets J- 1 C J- there is an FPRAS for Holant(J r ') 

The reason to take T' to be finite is to make sense of the computational problem 
Holant(J r '). As in Theorem [TJ if one is careful about how the input is specified, it is 
also possible to allow infinite T' in some cases. 



1.3 Matching circuits 

In Section [7] we will consider a natural type of gadget for reducing Holant problems to 
^PerfMatch, the problem of counting the number of perfect matchings in a graph. 

#PerfMatch is, famously, #P-complete even when restricted to bipartite instances [2"5] . 
This suggests that there is no efficient exact algorithm, leaving the question of whether there 
is an approximation algorithm. A major result in this direction is that there is an FPRAS for 
#Perf Match restricted to bipartite graphs 20J . Our study of matching circuits is an attempt 
to identify which Holant problems reduce to ^PerfMatch in the sense of expressibility. 

Consider a clique of order four, where at the z'th vertex we attach an "outgoing" edge 
di. For each of the sixteen possible subsets M C {^1,^2,^3,^4} of the outgoing edges, 
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we can count the number F(M) of ways to add internal edges to F to obtain a perfect 
matching. Because the clique of order four has 3 perfect matchings, we have -F(0) = 3, 
while F({dx}) = and F{{d 1 ,d 2l d 3 ,d 4 }) = 1. 

We will say that a function F : {0, 1} J — > Q-° has a matchings circuit if there is 
a similar graph fragment, with outgoing edges J, and such that F(x) is the number of 
perfect matchings containing the outgoing edges {i € J \ xi — 1}. Substituting each 
vertex of a Holant({F}) instance by the graph fragment gives a reduction from Holant({F}) 
to #PerfMatch. Actually, in Section [7J following Jerrum and Sinclair we will allow non- 
negative edge-weights and a "fugacity" at each vertex, because these do not add any more 
computational power; the important property is: 

Proposition 5. If T is a finite set of weight-functions that have matchings circuits, then 
Holant (J") < A p #Perf Match. 

Here <ap denotes a type of approximation-preserving reduction used to study the 
relative complexity of approximate counting problems - see Section 12.21 In particular, 
if Holant(J") <ap #PerfMatch and #PerfMatch has an FPRAS then Holant (J") has an 
FPRAS. The main result is the following theorem. 

Theorem 6. Let F : {0, l} 3 — > Q-°. The following are equivalent: 

1. F is windable 

2. For all x\, X2, x 3 € {0, 1} we have 



3. F has a matchings circuit 

Theorem[6]gives a class of problems that reduce to counting perfect matchings. For exam- 
ple, the Holant problem allowing only the relation {(0, 0, 0), (1,0,0), (0, 1,0), (1,0, 1), (0, 1, 1)} 
reduces to #PerfMatch, but is not known to have an FPRAS. 

1.4 Related work 

A matroid is sbo (strongly basis orderable) |llj if for all bases A and B there is a bijection 
7r: A\B -> B\A such that for all X C A \ B the set (A Un(X)) \ X is a basis. Bouchet and 
Cunningham generalised the sbo property as linkability for the class of even delta-matroids, 
and showed that this class is closed under an analogue of circuits [SJ. These conditions are 
just windability over the two-element Boolean semiring (B = {0, l},max, min), for the set 
of bases when considered as a function {0, 1} J — > B, by taking the characteristic vector of 
characteristic vectors of bases. Gambin used the sbo property to approximately count the 
number of bases in certain matroids 15J. 

Valiant |26j introduced matchgates and matchcircuits, which are similar to matchings 
circuits but give efficient exact algorithms. Matchcircuits can be understood as planar 
graphs with edge-weights, with no restriction to non-negative numbers. Cai and Choudhary 
characterised the expressibility of matchgates 0. The name "matchings circuits" used in 
this paper is meant to suggest a version of matchgates. 

The focus on (the negative side of) expressibility for approximate counting problems ap- 
pears in [S] , where logsupermodular functions are shown not to express non-logsupermodular 
functions in the context of #CSPs. 

Yamakami |2Z] and the current author [53] have given partial classifications for classes 
of Holant problems. The bulk of these results deal with intractability: reductions from a 
named problem such as #SAT to a given Holant problem. The focus of the current paper is 
on tractability: either in the absolute sense of an FPRAS, or by reductions to #PerfMatch. 



+ 



+ 



< 



F(x 1 ,x 2 ,x 3 )F(l - xx, 1 - x 2 , 1 - x 3 ) 
F(x\,x 2 , 1 - x 3 )F(l -xi,l- x 2 ,x 3 ) 
F(xi, 1 - x 2 ,x 3 )F(l - xi,x 2 , 1 - x 3 ) 
F(xi, l-x 2 ,l- x 3 )F(l - xi,x 2 ,x 3 ) 
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A related #CSP with mixed signs appears in the context of the Tutte polynomial. By the 
proof of |16[ Lemma 7], the following problems are equivalent in the sense of approximate 
counting, for any fixed y < — 1. 

• #Perf Match 

• #CSP({FJ) where B y : {0, l} 2 -> is defined by B y (0,0) = B y (l,l) = y and 
S y (0,l)=fl tf (l,0) = l 

• evaluating the Tutte polynomial at the point {x,y) where (x — l)(y — 1) = 2 

1.5 Outline 

In Section[3J we adapt the conductance argument of Jerrum and Sinclair to "even-windable" 
functions, which are a slightly simpler version of windable functions. We study windable 
functions in Section UJ We study strictly terraced functions in Section [5j In Section [5] we 
establish Theorem [1] and Theorem UJ Finally, in Section [7] we discuss matchings circuits and 
establish Proposition [5] and Theorem [6] 



2 Preliminaries 

A configuration of a finite "indexing" set J is a function x £ {0, 1} 7 . A weight- function 

is a function F : {0, 1} J -> A copy of F is a function G : {0, I} 7 -> of the 

form G(x) = F(x o 7r) for some bijection it: J — > I. We will use bold face to distinguish 
between sets S C J and the characteristic vector S. Similarly the bold version of a relation 
R C {0, l}" 7 is the corresponding zero-one-valued weight-function. 

We will not distinguish between {0, lj-'t 1 '— > fe } and {0, l} fc , or between {0, l} 1 and {0, 1}. 
Also, we will sometimes allow indexing sets to be partially enumerated in a certain way. 
This is for notational power: the enumerated indices are easy to refer to explicitly, while the 
unenumerated indices are easy to fix. For all positive integers k and all finite sets J, when 
k + J is used as an indexing set it means the disjoint union of {I, . . . , k} and J. Elements 
of {0, l} fe + J will be denoted by {x\, . . . , Xk] y) where x\,. . . ,Xk £ {0, 1} and y £ {0, 1} J . 

For all sets I C J, all configurations p of I and all configurations x of J \ J, let (x, p) 
denote the unique common extension of x and p to a configuration of J. The pinning 
of a weight-function F : {0, l} 7 Q-° by p £ {0, 1} 7 (I C J) is the weight-function 
F' : {0, 1} J \ 7 -> defined by F'(x) = F(x,p). 

The distance \{i £ J \ Xi ^ between two configurations x,y £ {0,I}' 7 will be 
denoted d(x, y). We say F is evenQif d(x,y) is even for all x, y with F(x), F(y) > 0. Define 
x © y £ {0, 1} J by (x y) 4 = + j/i (mod 2). For all x £ {0, 1} 7 define x £ {0, 1} J by 
^ = l-x h and for all F : {0, 1} J define FF : {0, 1} J ->■ by FF(x) = F(x)F(x). 

For all F : {0, 1} 7 Q-° and y £ {0, 1} 7 define the flip of F by y to be the weight- 
function F' : {0, 1} J defined by F'(x © y) = F(x) for all x e {0, 1} J . For alii € J 
define £ {0, l}' 7 (where J is implicit) to be the characteristic vector of {i}. 

For all finite sets J define 

Evenj = {x £ {0, l}' 7 | a;, is even} 
Odd,/ = {x G {0, l} 7 | x i is odd i 
NAEj = {x G {0, 1} J | 1 < ^ < | J\ - 1} 

EvenNAEj = Evenj n NAEj 

Even j and Oddj are parity relations. The last relation EvenNAEj is only used for 
calculations (and only with \J\ even). 

1 This may be confusing terminology - an even function can be non-zero on vectors of odd Hamming wei 
But it is a common definition for delta-matroids. 
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Figure 2: Terminology for graph fragments. 



2.1 Circuits 

In this paper, circuits are a type of graph equipped with weight-functions at each vertex, 
and allowing external edges. A little care is needed to allow self-loops and asymmetric 
weight-functions. 

A graph fragment G is specified by: 

• a set J whose elements are called incidences 

• a set V G of vertices, and sets J G , v £ V G , that partition J G 

• a set A G C J G whose elements are called external edges 

• a partition E G of J G \ A G into pairs called internal edges 
See Figure^ 

A circuit <j> is graph fragment equipped with a constraint Ftf : {0, l} Jv — > Q-° for 
each vertex v. We can also use a relation R C {0, 1} J « as a constraint by taking Ff = R. 
We will drop the superscript <j) where there is only one graph or circuit in context. 

G is closed if it has no external edges. Standard graph-theoretic terminology extends 
to graph fragments. In particular we will refer to connected graph fragments. An edge is 
either an internal edge or an external edge. A vertex v and an internal edge e are incident 
if J v intersects e. If an internal edge e is uniquely identified by the vertices u, v it is incident 
to, we will denote e by uv. 

Given a circuit 0, for any configuration x of J, 

• x is an assignment (with respect to E) if xi = Xj for all {i, j} € E. 

• x.\j v denotes the restriction of x to J v . 

• The weight of x is wt^x) = H veV F v (x\ Jv ). 

Fhc weight- function of cf> is the function [0] : {0, 1} A defined by 



where the sum is over extensions of x to assignments x' : {0, 1} J — > Q-° with respect to E. 
If a weight-function F is equal to [0] , we will say that F has the circuit <fi. 

Another way to think of a circuit is as a "read-twice pps-formula" , a special case of the 
pps-formulas of :5j. For example, consider an equation 



[0](x)=^wt <A (x') 



(xe{0,l} A ) 



X' 



1 1 



y =o z =a 
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Note how on the right-hand-side, each bound (summed) variable appears exactly twice, and 
each free (unsummed) variable appears exactly once. Any equation of this form defines a cir- 
cuit in a natural way: incidences correspond to the variable occurrences x, y, y, z, z; vertices 
correspond to terms G(x,y),G(y, z), H(z); the sets J v are scopes for each term; external 
edges correspond to free variables; and internal edges correspond to summed variables. 

For any partition £ of a finite set J into pairs, for all non-negative integers k, a k- 
assignment with respect to E is a configuration x of J such that Xi = Xj for all but exactly 
k pairs {i,j} £ E. So an assignment is a O-assignment. For all closed circuits <f> and all 
integers k > define 

z h {<t>) = wt ^ x )- 

fc-assignmcnts x 

So Zo((f)) is just [0] (evaluated on the empty configuration). 
2.2 Computational definitions 

A counting problem is a function / taking instances (encoded as strings over a finite alpha- 
bet E) to non-negative reals. A randomised approximation scheme for / is a randomised 
algorithm that takes an instance x and error parameter e > and returns an approximation 
Z to /(x) satisfying 

Pr [e~ e f(x) <Z< e e f(x)] > 3/4. (2) 

A fully polynomial randomised approximation scheme (FPRAS) for / is a randomised 
approximation scheme that runs in polynomial time in \x\ and e . (To be concrete, we can 
require the error parameter to be specified by a binary integer e _1 .) 

Let / and g be counting problems. A randomised oracle algorithm A meeting the follow- 
ing conditions is an approximation-preserving reduction from / to g, and if such a reduction 
exists we write / <ap <?■ 

A takes inputs (w, e) where w is in the domain of /, and e > 0. The run-time of A is 
polynomial in \w\ and e _1 and the bit-size of the values returned by the oracle (this avoids 
requiring that the oracle gives concise responses). The oracle calls made by A are of the 
form (v, 6), where v is an instance of g and 6 > is an error parameter, such that \v\ + S^ 1 
is bounded by a polynomial in \w\ and e _1 (depending only on .4). If the oracle's outputs 
meet the specification of a randomised approximation scheme for g, then A is a randomised 
approximation scheme for /. 

The above definitions are based on [13] ; the main difference is that we allow non- integer- 
valued problems. 

For any finite set T of weight-functions define the following counting problem. 
Name Holant(J r ) 

Instance A closed circuit <j) using copies of weight-functions in T 
Output [</>] 

Since T is finite, it is not particularly important how the functions Ftf are specified. 
For concreteness: F$ should be specified by an index i into a fixed enumeration J- = 
{Fi, . . . , F\jr\}, along with a bijection from J$ to the indexing set U of Fi : {0, — > Q-°. 

By substituting circuits, if F has a circuit using copies of weight-functions from a finite 
set J 7 , then Holant^U {F}) <ap Holant^J"). This justifies the focus on expressibility in 
this paper. 

3 Even-windable functions 
3.1 Idea 

Windability is an abstraction of a property of the distribution of perfect matchings in a 
graph with external edges. We will illustrate the idea briefly by the arity 4 case, where 
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Figure 3: An example of constructing perfect matchings by symmetric differences. From left to 
right, M, M', MAM' (with P drawn in thick solid grey), MAP, and M'AP. 

windability is already used implicitly in |19j . But higher-arity conditions are important for 
showing that windability is preserved by circuits. 

Consider a graph G with four external edges ei,e2,e3,e4. For all X\, X2, x$, X4, let 
F{x\, X2, X3, X4) be the number of perfect matchings in G that include the outgoing edges 
{ei I Xi — 1}. So F(0, 0, 0, 0)F(1, 1, 1, 1 ) is the number of pairs of perfect matchings 
(M, M') such that M includes all the external edges and M' includes none. But for any 
such pair (M,M'), the symmetric difference MAM' consists of cycles and paths, and the 
path starting at e\ ends at either e%, 63, or e±, depending on the choice of [M, M'). Thus 
F(0, 0, 0, 0)F(1, 1, 1, 1) splits into three terms. Denote these by B((0, 0, 0, 0), (f , 1 , f , 1 ), M) 
where M is a partition of {1,2,3,4} into pairs: either {{1, 2}, {3, 4}} or {{1, 3}, {2, 4}} or 
{{1, 4}, {2, 3}}. We can similarly define B((l, 1, 0, 0), (0, 0, 1, 1), M) for example. 

When MAM' contains a path P from e± to e2, the sets MAP and M'AP are also 
perfect matchings - see Figure |3] The only external edges in MAP are and e±, while the 
only external edges in M'AP are e x and e 2 . Thus B((0, 0, 0, 0), (1, 1, 1, 1), {{1, 2}, {3, 4}}) 
equals B((l, 1, 0, 0), (0, 0, 1, 1), {{1, 2}, {3, 4}}. 

In this section, for simplicity, we will consider only even functions. 

3.2 Definition 

For any configuration x G {0, 1} J define Ai x to be the set of partitions of {i G J | Xi = 1} 
into pairs. In particular, if J2ipj x i IS °dd then A4 X = 0. 

A function F : {0, 1} J — > is even-windable (with witness B) if there exist values 
B(x,y,M) > for all x, y e {0,1} J and all M e M x&y , i.e. all partitions M of the set 
{i € J I Xi ^= yi} into pairs, satisfying: 

EW1. F(x)F(y) = E Me M x9y 7, M) for all xje {0, 1} J , and 

EW2. B(x, y, M) = B(x © S, y © S, M) for all x, y e {0, 1} J and all S G M G M x ® y . 

Note that in the second condition, 5 is a pair {i,j} in M: we are swapping the values of 
Xi and yi, and swapping the values of Xj and yj. By swapping a sequence of pairs, EW2 is 
equivalent to 

EW2'. B(x, y, M) = B(x © Si © • • • © S fc , y © Si © • • • © S fc , M) for all x, y G {0, 1} 3 and all 

Sx,...,S k e M G Mxffiy 

3.3 2-decompositions 

Using pinnings, the even-windability conditions can be stated in a form that is sometimes 
easier to check. A function H : {0, 1} J — > Q-° has a 2-decomposition if there are values 
D(x, M) > 0, where x ranges over {0, 1} J and M ranges over partitions of J into pairs, such 
that: 

1. H (x) = -D(x, M) for all x, where the sum is over partitions of J into pairs, and 

2. D(x, M) = D(x © S, M) for all x, M and all S £ M. 

In particular if | J| is odd then the first condition forces H to be identically zero. 

A function F is even-windable if and only if for all pinnings G of F the function GG has a 
2-decomposition. For the forwards direction, given a witness B that F is even-windable, for 
each / C J and each p G {0, 1} 1 define D p (x, M) = 5((x,p), (x,p),M) for all x G {0, l} 7 ^ 



9 



to obtain a 2-decomposition D p of the pinning of F by p. For the backwards direction, 
for each / C J and each p € {0, l} 1 , pick a 2-decomposition D p of the pinning of F by 
p. For all x, y € {0, 1} J , define -B(x, y, M) = D p (x', M) where p is the restriction of x to 
{i € J | Xi = jji} and x' is the restriction of x to {i G J \ Xi =/= Hi}- Then B witnesses that 
F is even-windable. 

Lemma 7. Let F : {0, 1} J — > Q-° with \ J\ < 3. If F is even then F is even-windable. 

Proof. Let G : {0, l} 1 be a pinning of F. 

If J = define -D(x, 0) = GG(x) where x S {0, l} is the empty configuration. Then 
GG(x) = J2m -C( x j ^0 where M ranges over the set {0} of partitions of I into pairs, so D 
is a 2-decomposition of GG. 

If |/| = 2, let z,j be the elements of / and define D(x., {{i, j}}) = GG(x). For all 
x e {0, l} 7 we have GG(x) = ^ M D(x, M) where M ranges over the set {{{«,j}}} of 
partitions of / into pairs, so D is a 2-decomposition of GG. 

If |/| is 1 or 3 then G(x) and G(x) cannot be simultaneously be non-zero because G is 
a pinning of the even function F, and J2iei Xi = 1^1 + Sie/(1 ~ x i) ( m °d 2). Thus GG 
is identically zero. There are also no partitions of / into pairs, so the empty function is a 
2-decomposition of GG. □ 

Lemma 8. Even,/ and Odd,/ have a 2-decomposition whenever \J\ is even. Evenj and 
Odd,/ are even-windable for any J. 



Proof. First consider Evenj. Fix a partition N of J into pairs. Define 
D(x, M) 



1 if M = N and J2ie j x i IS even 
otherwise. 



Then for all x € {0, 1} J we have Evenj(x) = ^ M D(x, M) (where the sum ranges over 
partitions M of J into pairs). Similarly for Oddj, define 



D(x, M) 



1 if M = N and j is odd 
otherwise. 



Then for all x e {0, 1} 3 we have Oddj(x) = J2m m )- _ 
Now consider a pinning G : {0, 1} K of Evenj or Oddj. If \K | is odd_thcn GG 

is identically zero, by the same argument used in Lemma [7] Otherwise, G = GG is cither 
Everts or Odd^, which we showed have 2-decompositions. □ 

The following argument gives a more difficult example of a 2-decomposition. It will be 
used later (in the proof of Lemma IT7|) to show that NAEj is windable. 

Lemma 9. Let J be an finite set with \ J\ even. Then EvenNAEj has a 2-decomposition. 
Proof. For each subset / C J of even order fix a partition Mi of / into pairs. Set 



L>(x, M) = 2 



-fe+2 



{/ C J I |/| is even, Xi and are odd, and M = Mj U M 

iGl i£.J\I 



S G Mj U Mj\j implies S Q I or S Q J\I. The conditions that J^iei ^ an< ^ Sie J\/ • T * 
are odd are therefore not affected by changing xtoxffiS. Thus D(x © S, M) = D(x, M) 
for all SeM. 

For any x, if EvenNAEj (x) = then D(x,M) = 0. If EvenNAEj (x) = 1, pick i, j 
with Xi = and Xj = 1. For each of the 2 fe ~ 2 subsets /' C J \ {i,j} there is a unique set 
I" ^ {* j j} such that the order of / = I' U J" is even and such that J2i<=i x i an d x i 
are odd. There are thus 2 k ~ 2 such subsets / for each fixed x, which gives J2m -^( x j M) = 1. 
So D is a 2-decomposition of EvenNAEj. □ 
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Figure 4: A circuit <f> and a link graph Le{M) for some M £ .A/f x where x is an assignment 
of E®. (In particular, the M drawn is a union of partitions M v £ A^ x |j •) Circles represent 
vertices of the circuit. Squares are incidences i £ of the circuit, and are filled black where 
Xi = 1. Elements of M are drawn as thick black lines. Elements {i, j} £ E are drawn as thin 
lines. 

3.4 Expressibility 

We will show that the weight-function of any circuit using even-windable functions is even- 
windable. We will use a certain graph associated with a choice of matching of incidences. 

Let M and E each be a set of disjoint pairs of some set. Define the link graph Le(M) 
to be the multigraph on the vertex set UseM & with edge set the disjoint union of M and 
{{i,j} £ E | i,j £ UseM (so edges in M C\ E give pairs of parallel edges in Lg(M)). 

Note that for each vertex i of Le(M), the degree of i is two if {i,j} £ E for some 
j £ UseM ^> an d otherwise i has degree one. So Le{M) consists of paths and cycles. 

We will use this graph later for the analysis of the near-assignments chain. For now, 
consider an assignment x of some circuit with internal edges E and external edges A, and 
let M £ Ai x . For any i not in A with Xi — 1, the unique j with {i,j} £ E satisfies Xj = 1. 
This means that i £ Uspm & nas degree 1 in Le(M) if and only if i £ A. So every path 
component of Le(M) ends in {i £ A \ Xi = 1}, and every such i is at the end of a path. See 
Figure |U 

Lemma 10. Let tf> be a circuit using only weight-functions that are even-windable. The 
weight-function of (j) is even-windable. 

Proof. Recall that V, J, J V ,A, E, F v denote vertices, incidences, vertices' incidences, external 
edges, internal edges, and vertices' weight-functions. For each v £ V pick a function B v 
witnessing that F v is even-windable. 

Consider a set M 1 of disjoint pairs of J. We will say that M' induces the set of pairs 
{hj} Q A such that there is a path from i to j in Le(M'). 

For all x, y £ {0, 1} A and all M £ M xS)y define 

B(x,y,M) = J2 E l[B v (x'\j v ,y'\j v ,M v ) 

x '>y' {M„ } inducing M vEV 

where: 

• J2 X > y < denotes the sum over assignments x' and y' extending x and y respectively. 

• S{Af„} inducing m denotes the sum over all choices of M v £ M[ X ' my >)\j v for each v £V, 
such that \J v€V M v induces M. 
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For all x, y £ {0, 1} A we have 



EE E l[B v (x'\ Jv ,y'\j v ,M v ) 

MSA^xffly x',y' {M„} inducing M vGV 

J2 E n^( x/ k>y'U' M ") 

x',y' {A/„} t>GV 

E Il F vWjMy , \j.) 

[0] (x) 14>] (y). 

Here £{m } denotes the sum over all choices of M v £ M-{m'®y')\ Jv for each v £ V: the sum 
over M eliminates the condition that U«gv M v induces M. 

Now fix x, y £ {0, 1} A and S — {i,j} £ M £ A4 x ©y For any choice of {M v } inducing 
M, there is a unique path component -P{m„} (also depending on x, y, S) from i to j in 
^s(Ut,ev Usgm By construction of the link graph, the vertices of P{m v } are a union 
of pairs S £ {J veV M v . In particular, for each v £ V, the intersection Pim v } H J„ is a union 
of pairs S £ M v . Using EW2' we have 

B(x,y,M) = ]T J2 Y[B v (x'U,y'U,M v ) 

x',y' {M„} inducing M veV 

= E E n^(( x '® p {^})i^'(y'® p {^})i.^' M -) 

x',y' {M„} inducing M veV 

= B(x©S,y©S,M). 

So B witnesses that [</>] is even-windable. □ 



}2 S(x,y,M) = 



3.5 The near-assignments Markov chain 

Throughout this subsection fix an even-windable weight-function F : {0, l}' 7 — > Q-° and a 
partition E oi J into pairs. This can be thought of as a circuit with one vertex. We will 
define and study the near-assignments Markov chain for (F,E). 

Set n — \J\. For each k > let fife denote the set of fc-assignments of J with respect to 
E that satisfy F(x) > 0. The state-space is f2 = f2o U 1^. The transitions are Metropolis 
updates to states at distance two. More specifically, the transition probability from x to y 
is defined to be 

r^min(l,F(y)/F(x)) ifd(x,y)=2 
P(x, y) = I 1 - £ £y : d (x,y<)=2 min(l, F(y')/F(x)) if y = x 
I otherwise. 

(We will not consider the initial state to be part of the Markov chain itself: the Markov 
chain is completely described by the matrix P £ M x .) Define a probability distribution 
7r on fi by 

tt(x) = F(x) / £ F(y) (x G fl). 
/ yen 

By abuse of notation we will also denote J2 xeX by 7r (^) f° r subsets X C f2. By 
adapting the arguments of |19j . we will show: 

Theorem 11. For all x. £ Q and all non-negative integers t, we have 



\ £ |P { (x,y) -Tr(y)] < ^(x)- 1 /2 ex p(-t 7 r(f7 ) 2 /n 4 ^ 
yen 
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□ 

-□ 
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Figure 5: The link graph Le(M) for some 2-assignment-matching M. Squares are elements of 
J. Elements of M are drawn as thick black lines. Elements {i,j} G E are drawn as thin lines. 

Here P* denotes the i'th matrix power. The factor of I is convention: the left hand side 
is called total variation distance. 

We will use a congestion argument, with the following definitions. A flow-path is a 
directed path 7 in the transition grapfQ equipped with a weight wt(7), and also equipped 
with a label so that a set of paths can include the same path more than once. A flow T 
from X C f2 to Y C is a set of flow-paths which each start in X and end in Y, satisfying 

wt(7) = ir(x)ir{y) for all x € X and y EY. 

paths 7 £ r from a: to y 

The congestion of a flow T is defined to be 

p(r) = max - — V wt( 7 ). 

transitions fx, y) 7T X P(X,V) z — ' 

W 1 , '" 1 er with 7 9 (x,y) 

Be aware that this is not the same as the definition in : we are using a set of weighted 
paths rather than an assignment of weights to paths. When applying the results of [24] in 
the proof of Theorem [Til we will need to sum the total weight along each (unlabelled) path. 

In the following arguments we will often use fc-assignments (with respect to E) of the 
form x©y for x S and y € f2fc 2 . Note that we do not require F(x) > for fc-assignments 
x, though we do require P(x) > for x 6 fifc (fife is the set of "satisfying" fc-assignments). 
For any non-negative integer k, a fc-assignment-matching (with respect to E) is a set M 
of disjoint pairs of J such that exactly k edges {i, j} £ E have exactly one endpoint, i or j, 
in UseM In other words, the characteristic vector of Uspm S is a /c-assignment. 

Consider a fc-assignment-matching M. Recall the definition of the link graph Le{M) 
given in Section 13.41 which consists of cycles and paths. For any i with i S UseM &i ^ [ie 
unique j with {i,j} € E satisfies j g Use a/ except for exactly k values i. Thus Le(M) 
has precisely k/2 path components. See Figure OH 

For all non-negative integers k define 



Z k = £ F(x). 



xef2 fc 

(This is Zk(4>) if we consider P as a one- vertex circuit 0.) 
Lemma 12. Z Z 4 < Z 2 Z 2 . 

Proof. We have 



z Q z, = F (x)F(y) = E E B ( x >y> M )- 

xGf2 x£f! MGM x ^y 

yGf2 4 ye0 4 



2 the directed graph with vertex set ft and an arc (x, y) whenever P(x,y) > 0. 
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For each 4-assignment- matching M, pick a path component of Le(M) and let Hm be the 
set of vertices of this component. Let B be a function witnessing that F is even-windablc. 
Each Hm is a union of pairs in M so by EW2', 

Z Z 4 =J2 Yl £(x©H M ,y®H M ,M). 

xGSlo MeM^y 

But (x H M , y © H M , M) determines (x, y, M), and x H M , y © H M e fi 2 . So 

x'eo 2 MeA4 x , ffiy , 
y'en 2 

x'eo 2 
y'efi 2 



□ 



Lemma 13. Assume Zq > 0. There is a flow Tq from Qo to ismj flow-paths of length 



at most n/2, and with congestion at most in 3 /7r(O ). 



Proof. We will first construct a "winding" enumeration S(M, 1), . . . , S(M, \M\) of each 0- 
or 2-assignmcnt-matching M. The property we will need is that for each < k < |M|, the 
characteristic vector of S(M, 1) U • • • U S(M, k) is a 0- or 2-assignment. 

First define the final pair T(M) — S(M, \M\) for all non-empty 0- or 2-assignment- 
matchings M as follows. If M is a non-empty O-assignment-matching (so Le(M) consists 
of cycles), pick any vertex i € Le{M). If M is a non-empty 2-assignment-matching, pick an 
cndpoint i of the unique path component in Le(M). In either case let j be the unique index 
with {i,j} e M, and set T(M) = {i,j}. In any case L E {M\{T(M)}) = L E (M)\{i,j} has 
at most one path component. 

So M\ {T(M)} is a 0- or 2-assignment-matching. By induction on \M\ — k define 

S(M, k) = T(M \ {S(M, k + 1), . . . , S(M, |M|)}). 

So M\{S(M, k+1), . . . , S(M, \M\)} is always a 0- or 2-assignment-matching. This completes 
the construction of S(M, k). 

Let B be a function witnessing that F is even-windable. Let Tq be the set consisting of a 
flow-path 7x,y,M for each x € fi and y £ Q and M £ M x( £ y , where 7 x , y ,M is the flow- path 

x = x © J M ,o -> x © J M ,i ► x © J M ,|M| = y 

equipped with weight B(x, y, M)/(Z + Z 2 ) 2 and label (x, y, M), where Jju,t denotes the 
characteristic vector of S(M, 1) U • • • U S(M, k). 

r is a flow from £1 to O because for all x e il and y e Q 2 we have 

S(x,y,M)/(Z + Z 2 ) 2 - F(x)F(y)/(Z + Z 2 ) 2 = 7r(x)7r(y). 
The congestion of Fo is 

p(r ) = max ; ; \-, wt(7) 

HV ' transitions (»,»') 7r(z)P(z,z') ^ , W7 

7 6 Tq with (z, z ) G 7 

But 7r(z)P(z, z') = ^min(7r(z),7r(z')), so 
p(r o ) < max — ^ wtf/y) 

z€fi 2 • 7T(Z) 

7 £ To with z e 7 



max ■ 



£ £ B(x,y,M) 



zef2 2F(z)(Z a + Z 2 ) 

yefi 2 w ith ze7x, y ,M 
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In the last summation, z 6 7 x ,y,M implies z = x © 3m,k for some fc, so by EW2' we have 
B(x, y, M) = B(z, z © w, M) where w = x © y. Thus, 



^ r °^ 1 z ^ 2F(z)(Z + Z 2 ) E EE 5(z,z©w,M) 

0- and 2-assignmcnts w x£f2 MEMv, 

with ze7 x ,x©w,M 

For each (z,w,M) with M e .M w , the only values of x such that z 6 7x,xew,Af are the 
|M | + 1 values z © J M ,o, . . . , z © Jm,|m| • Thus, 

n 2 

p(Q ) <max— — — — — (|M| + 1) £j S(z,zffiw,M) 

zef2 Al 1 (z)(Zn + Zo) * — ' * — ' 

wv ; 0- and 2-assignmcnts w MeM w 

n 2 (n/2 + 1) v-^ , , . 

= max ,\ \ > F(z)F(z©w) 

zen 2F(z) Z + Z 2 ) ^ V ' V 7 

0- and 2-assignmcnts w 

Using z © w € fi U ^2 U fi 4) 



2 Z + Z : 



2 



li Z2 — then by Lemma [T2l we also have Z4 = 0, so the congestion is at most n 3 /2. 
Otherwise by Lemma H"2l we have Z4/Z2 < Z2/Z0 and 

^0 © ^2 © ^4 ^4 Zo 1 Zr\ , , _ , 

u < 1 + — < 1 + — = 1/ — = 1/tt(O )- 

Zq + Z2 Z2 Zq Zq + Z2 

□ 

Lemma 14. Assume Zq > 0. There is a flow T from to Q, using flow-paths of length at 
most n, and with congestion at most n 3 /tt(Qq) 2 . 

Proof. As in |20j . we will randomly route through fio- 

For each pair of flow-paths g, g' starting at the same state y, construct a flow-path 7(3, g') 
by appending g 1 to the reverse of g, and assigning a weight of wt(<?) wt((/')/7r(y)7r(f2o), and 
assigning the label (g,g')- Let To be the set of flow-paths given by Lemma [T51 Let ro(y,x) 
denote the set of flow-paths in To starting at y and ending at x. Let T denote the set of 
flow-paths 7(3, g') with g € To(y, x) and g' £ To(y, z) for some x, z €1 f2 and y £ fly. Then 
r is a flow from ft to f2 because for all x, z G O we have 

]T ]T wt( 5 )wt( 5 ')My)7r(n )= 7r(x)7r(y)7r(y)7r( Z )/7r(y)7r(a ) 

yeO 9 6r (y,x) yeo 

s'er (y,z) 

= 7r(x)7r(z). 

Letting (w, w') denote an arbitrary transition, the congestion of T is 



p(T) = max — — : — ; 



wt(<?) wt(g') 

fw.w') tt(w)P(w, w') ^— ' 7r(y)7r(f2n) 
v ; v y v ' ; x.zen ger ( y> x) yy ' y 01 

yen s ' e r (y,z) 

such that (w, w') £ 7(9,9 ) 



15 



By symmetry, 



p(Y) = 2 max ; ; V V 



E E 



wt(<?) wt(g') 

fw.w') tt(w)P(w, w') ^— ' ^— ' 7r(y)7r(On) 
yen 9 ' £ r„(y,z) 

such that (w, w 7 ) £ 



(w,w') tt(w)P(w,w') ^ 7r(y)7r(r2 ) 

x.ztii g £lo(y,z) 

yeO such that (w, w') e g' 

= 2 max - — V V wt( 5 ')M^o) 

(w,w') tt(w)P(w, w') ^ ^ w >i \ u> 

yeS '°such that (w,w') e g' 

= 2p(T )/7r(fio) 
< n 3 /7r(^o) 2 

by Lemma [T51 □ 

The remaining task is to relate the congestion to Markov chain mixing. 
Theorem llll For all x £ f2 and a/Z non-negative integers t, we have 

\ £ |P 4 (^y)-T(y)l < j7r(x)- 1 /2 exp (-t 7r (Q )2/n 4 > 
yen 



Proof. The transition matrix P is reversible relative to 7r: it obeys the detailed balance 
condition 

7r(y)P(y, z) = tt(z)P(z, y) for all y,ze!l. 

We have 

P(x, x) > 1 - 4- f ™ I > V" for all xe!!. (3) 
n 2 \2 J 

In particular, the Markov chain is aperiodic. Also, by Lemma IT4l there exists a flow Y from 
fl to f2, which implies that that the Markov chain is connected. This allows us to use the 
results from [23] and [T^]. P has eigenvalues 

f = A > Ai > ...A|n|_i > -1- 

By setting 7(7) in [MJ Corollary 6'] to be the sum of the weight of flow-paths in Y whose 
underlying directed path is 7, we have 

Ai < 1 ' 



p(Y)n 
< 1 - 7r(fi ) 2 /" 4 - 

using Lemma [T31 By (J3|) and equation 1 of [T7] we have 

-A| |-i <l-2/n< A x . 

By [TH Proposition 3], 



i^|P*(x,y)-7r(y)| < \ J^—^ max(A ls -A|n M ) 
yen V ^ ' 

< ivr(x)- 1 / 2 exp(-^(r!o) 2 /« 4 )- 



□ 
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4 Windable functions 



In this section we extend the analysis of even- windable functions to windable functions. The 
definition of windability is a natural extension of even-windability, but turns out not to give 
much extra generality. 

For all F : {0, 1} J -> define F e : {0, 1} 1+J by 

i^;x) = \y + iseven (p € {0,1},* €{0,1}'). 

I otherwise 

Lemma 15. F : {0, 1}' 7 — > Q— is windable if and only if F® is even-windable. 

Proof. (=>) Pick an ordering of J. Consider a partition M of a subset I C J into singletons 
{<zi}, . . . , {dfc} and pairs Si,..., Si- Define n{M), when |/| is even, to be the union of 
{Si,..., Si} with a partition (depending only on M) of {ai,...,afe} into pairs. Define 
(i(M), when |/| is odd, to be the union of {Si, . . . , Se} with a partition of {I, ai, . . . , a^-} 
into pairs. Let B be a witness that F is windable. For all (p; x), (g; y) € {0, 1} 1+J and all 
(p;x)©(q;y)7 define 

{J2 B(x,y,M') if p+J2i & j x i and 9 + X^ e j2/* are even 
M' :/1 (Af')=M 
otherwise 

For all S€M = fi(M'), if we let S" = S\{1} then S(x ® S', y © S', M') = B(x,y,M'). So 
i?' witnesses that F$ is even-windable. 

For all sets M of disjoint pairs of 1 + J define u(M) to be {S \ {1} | 5 G M}. Let 
i? be a witness that i 7 ® is windable. For all x, y g {0, 1} J and all M £ Ai' xS)y , define 

l 

B'(x,y,M)= Yl B((p; X ),(q;y),M') 

p,q=0 M' ■ v(M')=M 

For all S e M = u(M'), let S' = S if |5| = 2 and S" = S U {1} otherwise. Then 
£(x © S', y © S', M') = B(x, y, M'). So 5' witnesses that F is windable. □ 

Lemma 16. Let <fi be a circuit using only weight-functions that are windable. The weight- 
function of <p is windable. 

Proof. Replace each constraint F v by (F v )^, rename the new incidences p v , v £ V{<p), and 
add a constraint Eveni + p where P = {p v \ v £ V{(j))}. This produces a circuit </> ffi with 
[</>©] = [4>}q. By Lemmas ITSlfTUl and E] we find that {(j)®} is even-windable. So by Lemma US] 
again, [</>] is windable. □ 

Lemma 17. For any J, the weight- functions Even,/, Odd,/, and NAE,/ are windable. 

Proof. By Lemma [5] there is a witness B that Even/ is even-windable. Extending B by 
setting 5(x, y, M) = for all M £ M' xey \ -M x © y we get a witness B' that Even/ is 
windable. Similarly, Odd/ is even-windable by Lemma El and it is therefore windable. 

For NAEj, by Lemma [T5l it suffices to show that the weight-function (NAEj) ffl is even- 
windable. Let I C 1 + J, let p £ {0, 1} 1 , let K = (1 + jy\ J and let G : {0, 1} K ->■ <Q^° 
be the pinning of (NAEj)^ by p. We wish to show that GG has a 2-decomposition. If \K\ 
is odd then GG is identically zero so has a 2-decomposition. We can therefore assume that 
| if | is even. 

Let c £ {0, 1} be equal to |J| modulo 2. (NAEj) e is the weight-function corresponding 
to the relation Eveni + ,/ \ X where X consists of two configurations at distance \J\ + c. 
Specifically, X contains the all-zeros configurations of 1 + J, and also contains (c; x) where x 
is the all-ones configuration. We first argue that in all cases, GG is either Even/c or Odd/f , 
or a flip of EvenNAE^. 
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Figure 6: A circuit with weight-function F with F(0,0) = 2 and F(l,l) = 1 and -F(0, 1) 
F(1,0) =0. Vertices represent "exact-one" constraints {(1, 0, 0), (0, 1, 0), (0, 0, 1)}. 



If YlieiPi i s °dd, then G takes the value 1 precisely on Oddx \ X' where X' consists of 
at most one configuration x £ Odd^. If X' = then GG = Odd^. If X' = {x} then GG 
is the flip of EvenNAE^ by x. 

If J2iei Pi- 1S even ) then G takes the value 1 precisely on Even^- \ X' where X' consists 
of at most two configurations in Even^-. If \X'\ < 1 we are done by the same argument as 
the previous paragraph: GG is either Even^ or a flip of EvenNAE^. If \X'\ = 2 then 
X' consists of two configurations x, y with d(x, y) = \ J\ + c. But \ J\ + c < \K\ < \ J\ + 1, 
and \K\ and |J| + c are both even, so \K\ = \ J\ + c. Thus y = x, and again GG is a flip of 
EvenNAE A -. 

By Lemma[5]the weight-functions Even^- and Odd^- have 2-decompositions. So we only 
need to check the last case where GG is a flip, by z £ {0, 1} K say, of EvenNAE^. Let D 
be a 2-decomposition of EvenNAE^ given by Lemma [9j Define -D'(x, M) = £>(x © z, M) 
for all x G {0, 1} K and for all partitions M of K into pairs. For all x £ {0, 1} K we have 
GG(x) = EvenNAE^ (x © z) = J2m D ( :>1 © z ' M ) = Em ^'( x ' M )> wh ere M ranges over 
partitions of K into pairs. So £)' is a 2-decomposition of GG. □ 



5 Strictly terraced functions 
5.1 Idea 

To apply Theorem [11] to Holant problems, the challenge is to find a class of circuits for 
which the ratio of the weight of 2-assignments to the weight of 0-assignments is polynomially 
bounded in the size of the Holant instance. 

The weight of 2-assignments of a closed circuit can be written in terms of the weight- 
functions obtained by breaking two edges; the challenge then reduces to trying to find a 
bound on the ratios F(x) /F(y) between the values in these weight-functions, when F(y) ^ 0. 

It is instructive to consider multiplication of two-by-two matrices. To see the relationship 
between multiplication of matrices and circuits (in the form of read-twice pps- formulas) , for 
matrices M with rows and columns indexed by {0,1}, define Fm ■ {0, l} 2 — > Q-° by 
F M(i,j) = M itj ; then F MN (i, k) = £\ F M (i,j)F N (j, k). 

Matrix multiplication can produce exponentially-large ratios: for any x, y > 0, we have 

fx n I/O™" 1 + • • • + 1)\ 
[0 1 J 

and x n /l is exponentially large if x > 1. 

In fact, the matrix ( q ? ) corresponds to the circuit depicted in Figure fusing "exact-one" 
constraints {(1,0,0), (0, 1,0), (0,0, 1)}, which can be used to construct counterexamples to 
the bound ((T|) on nearly perfect matchings 0. 

We might guess that exponentially-large ratios can only be produced by matrix multi- 
plication when the zero entry in the matrix is surrounded by values that are different. And 
indeed this property of being "strictly terraced" turns out to give some control over ratios. 
For strictly terraced functions, the worst ratio in a weight-function is bounded by the sum 
of the worst ratios that can be obtained by mixing the individual functions with parity 
relations. 



fx y\ 
U 1 
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5.2 Definitions 



A function F : {0, 1} J -> Q-° is strictly terraced if 

F(x) = => F(x ej) = F(x © e^) for all x e {0, 1} J and all i, j e J. 

For all weight-functions F that are not identically zero, a parity-weight-function of 

F is a constant multiple of the weight-function of a circuit using one F constraint and such 
that all other constraints are parity relations. Define 

F'\ {0, 1} -» Q-° is a parity-weight-function 1 
of F with > J ' 

We extend 9 to all weight-functions F by setting 9(F) — if F is identically zero. 

We can show that 9 is well-defined using the following operation. For any circuit <f> and 
any internal edge e £ between incidences i u € J u and i v £ J v , with u,v £ V (not 
necessarily distinct), define the contraction of by e to be the circuit <fi' obtained by 
replacing u and v by a vertex w with incidences ( J u U J v ) \ {i u , i v } and equipping w with 
the weight- function of the circuit with constraints F u and F v , edge e, and external edges 
(J u U J v ) \ \imi>v}- 

Lemma 18. Let <f> be a connected circuit whose constraints are all parity constraints. The 
weight-function of (ft is a constant multiple of a parity constraint. 

Proof. By induction on the number of edges of cf), it suffices to show that contracting a single 
edge of (j) leaves only parity constraints (up to multiplication by constants). 

Consider the case that the edge {i, j} goes between distinct vertices, which are equipped 
with a Even^jy/ constraint and a Evenj^yj constraint. Then contraction gives a copy of 
Even/yj, because Even/uj(x, y) = J2t Eveni + /(t; x)Eveni + j(r; y) for all x S {0, l} 1 and 
all y e {0, 1} J . Similarly Odd^jyj and Odd^juj produce Even/yj, while Even^jy/ and 
Oddjjjuj produce Odd/yj. 

If the edge {i, j} is a loop on a vertex with constraint Even^j jyj, contracting {i,j} 
produces 2Even,/. Similarly Even^ ^uj produces 20ddj. □ 

Note that contracting an edge does not affect the weight-function of a circuit. By con- 
tracting edges between parity relations, the circuits appearing in the definition of a parity- 
weight-function can be rewritten not to use any edges except external edges and edges 
incident to the F constraint. For fixed F there are therefore a finite number of equivalence 
classes of parity-weight-functions F' : {0, 1} — > Q-° with F'(l) > 0, under the equivalence 
relation of multiplication by constants. Thus the maximum in the definition of 9(F) is 
taken over a finite set, which can be seen to be non-empty if F is not identically zero (if 
F(x) > for some x e Even./ then the function F' : {0, l} 1 ->• Q-° defined by F'(t) = 
Exe{o.i}-' Oddi(t)F(x)Evenj(x) satisfies F'(l) > 0, and if F(x) > for some x e Oddj 
then the function F' : {0, l} 1 -> defined by F'(t) = Exe{o,i} J Oddi(t)F(x)Oddj(x) 
satisfies F'(l) > 0). 

Note that if G is a parity- weight-function of F, then 9(G) < 9(F). In particular if G is 
a pinning of F then 9(G) < 9(F). Also, since the disequality relation Odd2 = {(0, 1), (1,0)} 
is a parity relation, it is not important that we took F'(0)/F'(1) rather than F'(1)/F'(0) 
in the definition of 9. 

5.3 Examples 

A relation R C {0, 1} J is coindependent if for all x e {0, 1} J \ R we have x ffi e^ e R 
for all indices i. For example, the disequality relation {(0, 1), (1,0)} is coindependent. Any 
coindependent relation R gives an example R of a strictly terraced weight-function. 

Lemma 19. For all finite sets J, the functions Evenj, Oddj and NAEj are strictly 
terraced. Also, 6>(Evenj) = 6 (Oddj) = 0, and 6>(NAEj) < 3. 



9(F) = max 



no) 
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Proof. The first statement follows from the fact that the corresponding relations are coin- 
dependent. To show ^(Evenj) = 0(Oddj) = 0, note that by Lemma IT51 a parity-weight- 
function of a parity relation must be even. 

Now we will show that 0(NAE,/) < 3. Consider a connected circuit <j> with one external 
edge, such that <f> uses one NAEj constraint, and all other constraints are parity relations, 
with no internal edges between parity relations (this is without loss of generality, because 
we can contract any such edge). Assume that [0] (0) and [0] (1) are non-zero. We will show 
that |0](O) < 3 

We can write 

xENAEj 

where R is an affine subspace of GF(2) 1+J . Since [</>] (0) and [0] (1) are non-zero, the sets 
Rq = {x | (0;x) 6 R} and Ri = {x | (l;x) 6 R} are non-empty. Since R is an affine 
subspace, |i?o| = so 

M (0)<|i?o| = |i?i|<M(l) + 2<3M(l). 

□ 

5.4 Properties 

An important property we will use is that a strictly terraced function F is either identically 
zero or its support {x | F(x) > 0} is coindepcndcnt. (If F(x) = and ^(y) > for some 
y, pick such a y with d — d(x, y) minimal. If d > 1, there are distinct indices i,j such that 
Xi 7^ yi and Xj ^ yj, so F(y © e,-) = F(y © © e^) = by minimality of d(x, y), which 
means F is not strictly terraced: F(y © e*) = but F((y © e,-) © e,) ^ F((y © e,) © e 3 ).) 

The Cartesian product of coindependent relations is in general not coindependent, for 
example {(0, 1), (1, 0)} x {(0, 1), (1, 0)} is not coindependent (set x = (0, 0, 0, 0) and i = 1). 
Thus the class of strictly terraced functions is not closed under taking weight-functions of 
disconnected circuits. 

Lemma 20. Let (f> be a connected circuit using strictly terraced weight-functions. Then [0] 
is strictly terraced. 

Proof. We will argue by induction on the number of internal edges of 4>. If there are no 
internal edges, then (j) consists of a single constraint using a strictly terraced function F, and 
[0] = F. Otherwise, pick an internal edge e. We wish to argue that the function created by 
contracting e is strictly terraced. There are two cases. 

(i) e is loop on a vertex v. 

Let F : {0, l} 2+,/ — > Q-° be a copy of F v , indexed so that the ends of e become 
enumerated indices. We wish to show that the function H : {0, 1} J — > Q-° defined by 

l 

tf(x)=]TF(M;x) (xe{0,l} J ) 
t=o 

is strictly terraced. Consider x e {0, 1} J satisfying H(x) — and let i, j G J. Since F is 
strictly terraced and F(0, 0;x) = F(l,l;x) = 0, we have F(0,0;x® e,-) = F(0,0;xffiej) 
and F(l, l;x © ej) = F(l, l,x © ej). Hence 7?(xffi e,) = H(x © ej). 

(ii) e is incident to distinct vertices u and v. 

Let F : {0, -> and G : {0, -> be copies of F u and F v respectively, 
reindexed so that the ends of e become the enumerated indices (and with / and J disjoint). 
We wish to show that the function H : {0, 1} /UJ -> Q-° defined by 

l 

H( X , y)=E F(t; x)G(t; y) (x G {0, l} 1 , y e {0, 1} J ) 

t=0 
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is strictly terraced. If F or G is identically zero then H is identically zero and therefore 
strictly terraced. 

Otherwise consider x £ {0, l} 7 and y £ {0, 1} J satisfying if(x, y) = 0. Since F and G 
have coindependent support and F(0; x)G(0; y) + F(l; x)G(l; y) = 0, there exists f € {0, 1} 
such that F(t; x) = G(l - <; y) = and F(l - i; x), y) > 0. For alH € / we have 

H{x © e*, y) = F(t; x © ei )G(t; y) = F(l - t; x)G(t; y). 
Similarly for i e J we have 

JT(x, y © a) = F(l - t; x)G(l - t; y © e 4 ) = F(l - t; x)G(i; y). 
Therefore for all z, j £ / U J we have if((x, y) © = jff((x, y) © ej). □ 

The following calculations bound ratios produced by certain circuits. 
Lemma 21. LetF: {0, 1} 1+J and G : {0, 1} J -> Define H(Q),H(1) by 

H{t)= ^(*;x)G(x). 

x£{0,l} J 

^4s5wme £/ia£ i* 1 and G are strictly terraced and that H(l) > 0. TTien 

H(Q)<WF) + ^(G))JJ(1). 

Proof. We will use induction on |J|. For the base case J = we have -ff(O) < 9(F)H(1) by 
definition of 0(F). So assume that J is non-empty. 

For each i £ J and each c € {0, 1} define Fi, c to be the pinning of F by taking i to c, 
and similarly define Gi jC to be the pinning of g by taking i to c, and define 

ffi,c(*)= J2 F,,il.xiG,,ix: (te{0,l}). 

xe{o,i} J \{*} 

Since pinnings are parity-weight-functions, 9(Fi }C ) < 9(F) and 9(Gi :C ) < 9(G). If there 
exists i £ J such that Hi q(1) and Hi t i(l) are non-zero, then by the induction hypothesis 
we have 

H(0) = H it0 (0) + Hi,i(fi) < (9(F) + 9(G))(H hQ (l) + H hl (l)) = (9(F) + 9(G))H(1). 

Taking a choice for each i, we may assume that there exists y £ {0, 1} J such that for all 
i £ J we have Hi^- yi (l) = 0. So for each i £ I the sets R — {x | F^x-^ljx) > 0} and 
S = {x | G^i-y^x) > 0} are disjoint. R and S are pinnings of coindependent relations, 
so they are coindependent. For all x £ R we have x ^ S, so x © e, € S for any i, and 
x © ^ R. Repeatingly this, we find that R consists of the configurations at even distance 
from x, and S consists of the configurations at odd distance from x. In other words, for 
each i £ J there exists Ci £ {0, 1} such that 

F(l, x) > <^=> Ci + Xj is even, and 

^ (xe{0,i} J ,^y,) (4) 

G(x) > ^=> Ci + ^ Xj is odd. 

For any i, j £ J there is some choice of x <= {0, 1} J with a;, ^ yi and Xj 7^ yj, so Cj = Cj. 
Thus there is a single choice of c such that (U]) holds for all i taking Cj = c: 



x) > c + xj is even, and 



G(x) > c + V" is odd 



(xe{0,i} J \W) (5) 
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For any x £ {0, 1} J with c + 2~^ig j x i even i an d anv distinct i, j £ J, we have F(l; x) = 
F(l;x (B ei (B ej) because F is strictly terraced and either F(l;x ® e^) or F(l;x ffi e.,-) is 
zero. Similarly for any x £ {0, 1} J with c + 2~Ziej odd, and any distinct i, j £ J, we have 
G(x) = G(x © ei © ej). This implies that there are constants A, fj, > such that 



(xe{0,l} J ). (6) 



If c + J2ie.J Vi ^ °dd then by (J6j) and ([5]), G is /iEvenj (if c = 1) or /iOddj (if c = 0). 
So G is a constant multiple of a parity relation. Considering if as the weight-function of 
a circuit with constraints F and G, we get H(0) < d(F)H(l) by definition of 9. We may 
therefore assume that c -I- J ^» ^ s even - 

Define F'(0), F'(l), G'(0), G'(l) by 

F'(i) = F(i;y) 

G'(t)= ]T Odd 2+J (t, c;x)G(x) 
xe{o,i} J 



F(l; x) = A if c + x i is even, and 
G(x) = /i if c + is odd. 



iG J 



so 



F'(0) = F(0;y) 
F'(l) = F(l;y)=A 

G'(0) = m2 |j|_1 
G'(l) = G(y) 

Since F' is a parity-weight-function of F we have F(0;y)/A < 6(F). Since G' is a parity- 
weight-function over G we have /i2l J ' _1 /G(y) < 9(G). For all x e {0, 1} J with c + x i 
odd, we have F(l;x) = and F(l;x© ei) = A (for any i £ J) and therefore F(0;x) = A 
because F is strictly terraced. So 

H(0) _ F(0,y)G(y)+A / ,2l^l- 1 <f)(p) , 
W) ~ AG(y) * + 9{G) - 

□ 



Lemma 22. Let (f> be a circuit using strictly terraced weight- functions. Then 

em) < J2 e ( p S)- 

Proof. We will argue by induction on the number k of constraints that are not parity rela- 
tions. The cases k = and k = 1 follow from the definition of 9. Components of a circuit not 
connected to the external edges simply contribute a constant factor to the weight-function. 
So for a given fc, it suffices to show that [0(0)] / [</>(l)| < X^ev* &(Fv) whenever: 

• <j) is a connected circuit with one external edge, with [0(1)] > 0, and 

• <j) uses strictly terraced weight-functions, at most k of which are not parity relations. 

For the k = 2 case, if there is a loop on a vertex v, contract it. This changes the 
weight-function F v , but the resulting weight- function is a parity- weight- function of F v , so 
this process does not increase J2vev ^(F v ). And F v is still strictly terraced by Lemma |2T)1 
Similarly, if there is an edge incident to distinct vertices u, v where F u is a parity constraint, 
contract that edge. The weight-function F introduced by the contraction is a parity-weight- 
function of F v , so this process does not increase Ylvev e(F v ). And again, F is strictly 
terraced by Lemma 1201 Repeating this process we end up with a circuit with at most two 
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Figure 7: An illustration of the correspondence between arbitrary configurations in a closed 
circuit and assignments in a modified circuit. Solid circles are arbitrary constraints, empty 
circles are copies of Even3. Thin black lines are incidences given the value 0, thick grey lines 
are incidences given the value 1. 



vertices. If there is only one vertex we can appeal to the k < 1 case, and otherwise we are 
done by Lemma [5TJ 

For k > 2, contract any internal edge. From the k < 2 case we know that J^vev @(Fv) 
has not increased. This process does not change the weight-function of <fi, and by Lemma 
[20] the constraint function introduced by the contraction is strictly terraced. □ 

Lemma 23. Let <j> be a closed circuit using strictly terraced constraints, and assume that 
Z ((j)) > 0. Then 



Proof. We will consider a circuit ip obtained by attaching Even3 to edges of <p as illustrated 
in Figure [7] In words: let J* be a disjoint copy of J, consisting of an clement i* for each 
i € J. Define ip to have incidences E U J U J* , external edges E, edges {i,i*} for each i 6 J, 
vertex set V U E, the same constraints at each v £ V, and Fuji = Even^* for all 

{i,j}eE. 

Zk((j)) is the sum of ftp] (x) over configurations x of E^ with X^ee-E* Xe = ^- By pinning, 
bounds the ratio F(x © e^/i^x) between the weights of neighbouring configurations of 
non-zero weight. Letting denote the all-zeros vector, for all i ^ j such that fipj (e^) ^ 0, 

If J?/'] (ej) = we have (e,: + a,-) = [^] (0) because fipj is strictly terraced. Thus 

Z 2 (4>) < Z Q (0)(\ E ^ max(l,0(V>)) 2 . 
The result follows by applying Lemmas and [22] □ 



6 Proofs of Theorem [T] and Theorem |4] 

Theorem [TJ #ParityNAE has an FPRAS. 

Proof. We are given a labelled graph, which is naturally a closed circuit <fi using constraints 
of the form Even,/, Odd,/, and NAEj. 

The decision problem, deciding whether Zo((f>) > 0, can be solved in polynomial time by 
Cornuejols' algorithm for the general factor problem [TO]. And degree-f parity relations can 
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be used to fix edges to take a particular value. This means that ^ParityNAE is self-reducible 
in the sense of [UJ Theorem 6.4]. So it suffices to give a fully polynomial almost uniform 
sampler (FPAUS): an algorithm that, when given an error parameter < 6 < 1 and an 
instance of ^ParityNAE corresponding to a closed circuit <f> with Zq{4>) > using parity and 
not-all-equal relations, outputs an assignment z satisfying 

i^|Pr(x = z)-wt (x)/Z o (0)|<£ 

X 

in time polynomial in the size of the input and in log<5 _1 . 

We will use the near-assignments chain to sample from assignments of <f>. Define F 1 : 
{0,1} J ^Q>° by 



f'(x)= n^(*D 



vev 

and 



F(x) = 



F '( x ) J2ieJ x * is even 
otherwise. 



By Lemma fTTl all the constraints of <fi are windable. By Lemma [TBI F' is windable. By 
Lemma [T5l is even- windable. But F is a pinning of (F')& (setting the parity bit to 

zero). A pinning of an even- windable function is even-windable - this is immediate from the 
characterisation in terms of 2-decompositions given in Section 13.31 

We will use the notation 7T, f2, Slo from Theorem II 11 for the near- assignment chain on the 
pair (F, E^). 

Recall from Lemma [TJ that 0(NAEj) < 3 and 6»(Even. 7 ) = 0(Oddj) = 0, and that 
all these weight-functions are strictly terraced. Let R = 3\V\ 2 \E\ 2 ; by Lemma |2"31 we have 
1/R < Zo(0)/Z 2 (0) < Z Q {<f>)/{Z {ct>) + Z 2 {<f>)) - tt(Qo). 

By Cornuejols' algorithm, mentioned above, we get an assignment y with [0] (y) > 
and in particular 7r(y) > 2 _ '- E L Applying Theorem 1111 by simulating the near-assignments 
Markov chain of (F, E) for t > (2\E\) A R 2 (log ^ + |£7| log 2) steps we can take a sample z 
from near-assignments of <f> such that 

^|Pr(x = z)-7r(x)|<5/2i? 

Thus 

i I Pr(x = z|z e fio) - F( X )/Z (cf>)\ < 6/2 

So we get an FPAUS by rejection sampling: run the simulation at least 2i?log| times 
and return the first sample in fio- The probability that this fails is small (at most (1 — 
i^(r2 )) 2i?1 °sl < 5/2). ' □ 

Theorem [4j Let T be the class of strictly terraced windable functions. Then 

• T is closed under taking weight-functions of connected circuits 

• T contains Everifc, Oddfc, and NAE^ for all k > 1 

• for all finite subsets J-' C T there is an FPRAS for Holant(J-"') 

Proof. The first statement is Lemma [TBI The second statement is Lemma 1 171 

For the third statement, given J 7 ', let T" = T 1 U {Eveni, Oddi}. We can use Eveni 
and Oddi to fix the value an edge takes, so Holant(J r ") is self-reducible. (There is a 
minor difference from [3JJ: we allow rational- valued functions. But this is not important.) 
For the decision problem we can use Feder's algorithm for coindependent relations [21 
Theorem 4]. Otherwise the argument proceeds as in the previous proof, taking R to be 
|£| 2 |V| 2 max/^jcv 0(F). We find that Holant(J r "), and therefore Holant(J r '), has an FPRAS. 

□ 
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7 Matchings circuits 

Define a matchings circuit G to be a graph fragment equipped with: 

• a non-negative rational edge-weight w(e) for each internal edge e 

• a non-negative rational fugacity X(y) for each vertex v 

Note that in this definition the external edges are not given weights. 

Let deg F (v) denote the number of edges in F incident to the vertex v. The weight of a 
set of edges F C A U E is: 



wt G (F) = 



if degp(v) > 2 for all vertices v 

XUd cgF ( v )=oMv)) (ILgWe)) otherwise. 



The weight-function of G is the function [G] : {0, 1} A — > Q-° where A is the set of 
external edges and 

[G] (x) = Yl wt ^ F )- 

FnA={e£A|a: e =l} 

As with circuits, if F = [G] we will say the F has the matchings circuit G. 
For all w > define Edge 10 : {0, l} 2 -> by 



Edge™(i,j) = 



1 

w 



where the matrix rows and columns are indexed by {0, 1}. For all A > and all finite sets 



J define Fugacity) : {0, 1} J -> by 



Fugacity^(x) = { 1 ifE ie j^ = l 
otherwise. 

Given G, define a circuit by equipping each vertex v with the function Fugacity then 

subdividing each edge e and equipping the new vertex with the function Edge 1 "^'. The 
circuit clearly has the same weight- function as the matchings circuit. So matchings circuits 
are just a special type of circuit. We will use the same notation and terminology. 



7.1 Example 

Proposition 24. For all finite sets J define OR,/ = {x e {0,f}" 7 | J2iej x i > 0}- Then 
ORj has a matchings circuit. 

Proof. We may assume J — {1, ... , k}. The matchings circuit is illustrated in Figure |S] 
Define F : {0, f } 3 -> by 

^,i,i)=(; ;).. 

(with rows and columns indexed from zero.) Each of the smaller boxes shown in Figure |S] 
have the weight- function F (where external edges are numbered from left to right). 
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Figure 8: 2 k 1 ORfc matchings circuit. Hollow circles are vertices with fugacity 1. All other 
vertices have fugacity 0, and all edges have edge- weight 1. 



For all xi, . . . , £ {0, 1}, 

[Gj(xi,...,x k )= ^2 F ( 1 ,Xi,yi)F(y 1 ,x 2 ,y2) ■■■F(yk-i,Xk,0) 




= 2 k - 1 OR k {x 1 ,...,x k ). 

So the weight-function of G is 2 fc ~ 1 OR/ t ,. To deal with the scalar multiple, add an 
isolated vertex with fugacity l/2 fe_1 . □ 

In particular, let k > I be odd. Then (OR&)© is a copy of EvenORfc + i where 
EvenORfe+i is defined as Even^+i n ORfe+i. By Lemma fT5l EvenORfc + i is even-windable. 
Thus EvenORfe + iEvenORfe + i = EvenNAE^+i has a 2-decomposition. This gives an 
alternate proof of Lemma |H] which, via Lemma H3 shows that NAEj is windable for all 
finite sets J. But this argument does not seem to show that NAEj has a matchings circuit. 



7.2 Approximate counting 

Define 

Name #PerfMatch 

Instance A simple graph G 

Output The number of perfect matchings in G 

The aim of this section is to establish Proposition[5j that Holant(J 7 ) <ap #PerfMatch for 
any finite set F of weight-functions of matchings circuits, showing that matchings circuits 
are a natural choice of circuit for ^PerfMatch. We will reduce via the following problem. 

Name #FugacityWeightedPM 

Instance A closed matchings circuit <f> where fugacities and edge-weights are given as ratios 

of non-negative integers specified in binary 
Output Zq{4>) 

The fugacities and edge- weights can both be simulated using matchings circuits. A 
similar reduction appears in |16j . 
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Figure 9: Gj^, with one path in the copy of GV 5 i labelled. All fugacities are 0, all edge- weights 
are 1. 

Lemma 25. There is a polynomial-time algorithm which, given non-negative integers p,q 
specified in binary, outputs a matchings circuit G Pi9 whose fugacities are all and whose 
edge-weights are all 1, and with two external edges such that 

{G P J{i,j)= (* °) for all i,j €{0,1}. 

where we consider the rows and columns of the matrix to be indexed from zero. 
Proof. See Figure [SI 

For all p > there is a unique binary expansion p = 2™ 1 + • • • + 2™ fc , with < n\ < ■ ■ ■ < 
rife. Define G p ,i in the following way Take two vertices s and t, each with one external 
edge. For each I < i < k, if n» = add an edge between s and t, and otherwise add a path 
between s and t of length 2ni — 1, which we can denote s — Vi,\, Vi,2, ■ ■ ■ , Vi^m = t, and add 
a parallel edge in the odd positions: between Uj,2j-i and Vi^j for each 1 < j < m. 

There is a unique perfect matching of G Pt \ that includes the external edges: it uses the 
edges in even position in each path, Vi^jVi^j+i for all 1 < i < k and all 1 < j < n,i. The 
perfect matchings of G v .\ that do not include the external edges are determined by a choice 
of 1 < i < k such that the i'th path uses edges in odd positions, and a choice of edge in 
each of the odd positions in this path; there are 2 Ui choices for each i. So G Pi i has the 
correct weight-function: [G p ,i] — (^J).^. 

Define H to be the circuit consisting of one vertex with fugacity zero, with two external 
edges. For q ^ 1 define G P:q to be serial composition of copies of G p .i, H, G qt %, and H, that 
is, we identify the second external edge of the i'th circuit with the first external edge of the 
(i + l)'th, for i = 1, 2, 3. The weight-function of G Ptq is then given by the matrix 

fp 0\ (0 1\ fq 0\ [0 1\ fp 
\0 l) \l Oj \0 l) \l Oj \0 q 

□ 

Lemma 26. Given a matchings circuit G for a weight-function F, we can efficiently con- 
struct a matchings circuit G' for (defined in Section^ in which every vertex has fugacity 
zero. Conversely, given a matchings circuit G for Fq, we can efficiently construct a match- 
ings circuit G' for F . 

Proof. For the first statement, pick an enumeration v\,...,v n of the vertices of G. Form G' 
as follows. For each 1 < i < n, add vertices dj, edges aibi,aiCi,biCi with edge-weight 

1, add an edge ViOi with weight X(vi), and if i < n add an edge Cibi+\ with edge-weight 
1. Set all the fugacities to zero and add an external edge at b\. See Figure ITD1 Consider a 
matching M C E of G. We will argue that there is a unique way to extend M to a perfect 
matching M' of G' . 

Let U = {i | dcg M (vi) = 0} be the indices of unmatched vertices. Let Mi = {a^ | i € 
U}. Note that the extension M' must include Mi, and if i ^ U we cannot have biCi £ M'. 
Consider the following subset P of external and internal edges: the external edge at bi, 
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Figure 10: Illustration of a matching circuit for built from a matchings circuit for F, as 
described in Lemma [26l 



edges biCi for all i <E U, and the edges biCii and a.;C; for all i £ U. So P is a path, except 
that at one endpoint, P has an external edge b\. Observe that there is a unique choice of 
perfect matching M 2 C P along this path: the odd-numbered edges starting from the end 
of P not incident to the external edge b\. (If P has an odd number of vertices then we get 
61 G M 2 , and otherwise b x g M 2 .) Define M' = M U M x U M 2 . Any extension of M to a 
perfect matching of G' would have to include M\, and hence M 2 , and so the extension is 
unique. 

This gives a weight-preserving bijection between matchings M of G and perfect matchings 
M' of G'. Since G' has an even number of vertices, the sets M' must include an even number 
of external edges. Thus \G'\ = P®. 

The converse is easy: given a matchings circuit G for P® , add a vertex of fugacity 1 to 
the external edge 1 to get a matchings circuit for P. □ 

Lemma 27. #FugacityWeightedPM < A p #PerfMatch 

Proof. Given a matchings circuit G\ with no external edges, we will construct a simple 
graph G with C \G\\ perfect matchings where C is an easily computed positive integer. 

By Lemma !2ol we get a matchings circuit G 2 such that \G 2 \ = [Gi]®- Deleting the 
external edge, we get a circuit G3 with [G3] = [GiJ. At each edge e of G3, we have integers 
Pe,q e such that the weight of e is p e /q e . Insert a copy of the circuit G Pi9 given by Lemma 

this produces a circuit G4 whose weight-function is C [G3] where C = n e 6_B G 3 9 e ' ano - 
where all fugacities are and all edge- weights are 1. Forgetting the fugacities and edge- 
weights we get a multigraph with G [G3] perfect matchings. To construct G, delete any 
loops and subdivide each edge into a path of length 3; this does not affect the number of 
perfect matchings. □ 

Proposition [5j If T is a finite set of weight-functions that have matchings circuits, then 
Holant(P) < A p #Perf Match. 

Proof. Pick a choice of matchings circuit Gf for each P G P. Given an instance of 
Holant(P), for each vertex v the function F v is a copy of some P G P; we can substitute 
Gf into tp at v. This process gives a matchings circuit G' with the same weight-function as 
ip. We can then appeal to Lemma B71 □ 

7.3 Expressive power 

Lemma 28. The weight-function of any matchings circuit is windable. 
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Figure 11: A weighted clique. All fugacities are zero, and w(uvi) = F(ei) for all i. The other 
edge-weights are to be determined. 

Proof. By Lemma [551 and Lemma [15] it suffices to show that the weight-function of any 
matchings circuit where all fugacities are zero is even-windable. For all w > 0, Lemma [7] 
implies that Edge™ is even-windable. For all A > and all finite sets J, consider a pinning 
G : {0, 1} 1 -> of Fugacity^. If (GG)(x) > for some x then £\ gJ x, and £ ieJ (l - Xi) 
are at most 1, so \I\ < 2. Thus GG has a 2-decomposition as in Lemma [7] □ 

To give circuits for low-arity functions we will apply linear programming duality in the 
form of Farkas' lemma. For a very short proof of Farkas' lemma, as well as a statement 
explicitly allowing a general ordered field, see pQ. For all x, G Q d denote the dot product 
Xi<j>i -\ h x d cj) d by x- (f>. 

Lemma 29. Let xi, . . . , X&, y S Q d . TTie following are equivalent: 

• y = cixi H h c fe x fe /or some c x , . . . , c k € Q-° 

• y • 4> > /or aZZ <f> 6 Q d £/ia£ satisfy Xi •</>,..., x& • > 

Lemma 30. Le£ F : {0, l} 4 -> Assume that F(et), F(eZ),F(e£),F(ei) are not all 

zero, and that F(x±, x 2 , ^3,^4) = whenever xi + x 2 + x 3 + x 4 is even, and that for all 
xi, X2, X3, X4 € {0, 1} we /iaue 

F(xi,a;2 ! a;3,a;4)F(l - xi, 1 - x 2 , 1 - x 3 , 1 - ar 4 ) 
< F(xi,x 2l 1 - 353, 1 - x 4 )F(l - xi, 1 - x 2 , 53 3 , ar 4 ) 

+ F(xi, 1 - x 2 , x 3 , 1 - x 4 )F(l - Xi,x 2 , 1 - 0:3,0:4) 

+ F(xi, 1-332,1- x 3 ,a;4)F(l - £1,2:2,2:3, 1 - x 4 ). 

Then F has a matchings circuit. 

Proof. We will construct values w{viVj) > satisfying 

F(ei) = F(e^)w(v 3 v 4 ) + F(e^)w(v4V 2 ) + F(e^)w(v 2 v 3 ) 
F(e 2 ) = F(e^)w(v4Vi) + F(e^)w(viv 3 ) + F(eT)w(v 3 v 4 ) 
F(e 3 ) = F(e^)w(viv 2 ) + F(W[)w(v 2 v 4: ) + F(e^)w(v4Vi) 
F(e 4 ) = F(eT)u!(v 2 v 3 ) + F(e^)w(v 3 vi) + F(e^)w(viv 2 ) 

(and also w{vtVj) — w(vjVi).) This suffices because then F = [G] where G is the weighted 
clique illustrated in Figure [TTJ We need to show that the vector 

y = (F( ei ),F(e 2 ),F(e 3 ),F(e 4 )) 
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is a non-negative linear combination of the vectors eiF(ej) + ejFfa) with 1 < i < j < 4. 
By Farkas' lemma (Lemma [23), it suffices to show that y • cj) > for all <p E Q 4 satisfying 

(/>iF(eJ) + <t)jF(e-) > for all 1 < i < j < 4. (8) 

If 4>i, (f>2, ^3, (/>i > we are done. Otherwise 0, < for some i. By assumption F(eJ) > 
for some j. If j ^ i, then (j5]) implies that tpjF^el) is non-zero. In any case -F(el) > 0. If 
i = l then 

F( ei )F(el) < F(e 2 )F(et) + F{e 3 )F(ej) + F(e 4 )F(ej) 
-«^iF(ei)F(er) < -<p 1 F(e 2 )F(e^) - c/> 1 F{e 3 )F(e^) - ^F(e 4 )F(el) 
-<^iF(ei)F(eT) < <f> 2 F(e 2 )F(et) + 3 F(e 3 )F(eT) + </> 4 F(e 4 )F(e7) 
-0i-F(ei) < 2 F(e 2 ) + <feF(e 3 ) + 4 F(e 4 ) 
Therefore y • </> > 0. By symmetry the other cases, i ^ 1, are similar. □ 

Theorem [BJ Le£ F : {0, l} 3 — > T/ie following are equivalent: 

1. F is windable 

2. For all xi,x 2 ,x 3 € {0, 1} we have 

F(xi,x 2 ,x 3 )F(l - xt,l- x 2 , 1 - x 3 ) 
< F(x 1 ,x 2 ,l-x 3 )F(l-x i ,l-x 2 ,x 3 ) 
+ F(xi,l-x 2 ,x 3 )F(l-xi,x 2 ,l-x 3 ) 
+ F(xi,l-x 2 ,l-x 3 )F(l-Xi,x 2 ,x 3 ) 

3. F has a matchings circuit 

Proof. For notational convenience, in the following argument we will use a particular copy 
of Fq. For all Xi,x 2 ,x 3 ,xn € {0,1} define 



F'{x-i,x 2 ,x 3 ,Xi) 



F(xi,x 2 ,x 3 ) if x\ + x 2 + x 3 + X4 is even 
otherwise. 



(1 => 2) Let B be a witness that F' is even- windable (using Lemma [T5|) . Let 
xi, x 2 , x 3 G {0, 1}. Let c € {0, 1} be the unique value such that x\ + x 2 + x 3 + c is even. 
Then (x\, x 2 , x 3 , c) © (1 — X\, 1 — £2, 1 — £3, 1 — c) = (1, 1, 1, 1). Note that 

*t( M ,i,i) = {{{L2},{3,4}},{{1,3},{2,4}},{{1,4},{2,3}}}. 

We have 

F(xi,a!2,a;3)F(l - x\,l - x 2 , 1 - x 3 ) 
= F'(xi 1 x 2 ,x 3 ,c)F'(l - xi,l - x 2 , 1 -0:3,1 - c) 
= B((x 1 ,x 2 ,x 3 ,c), (1 -xi,l -x 2 ,l -x 3 ,l - c), {{1,2}, {3, 4}}) 

+ B(( Xl ,x 2 ,x 3> c), (1 - x u 1 - x 2 , 1 - x 3 , 1 - c), {{1,3}, {2, 4}}) 

+ B((aji, x 2 ,x 3 , c), (l-x u l- x 2 , 1 - ar 3l 1 - c), {{1,4}, {2, 3}}) 
= B((xi,a;2, 1 - £3, 1 - c), (1 - xi, 1 - x 2 ,x 3 ,c), {{1,2}, {3,4}}) 

+ B{{xi, I - x 2 ,x 3 ,l - c), (1 - xi,x 2 ,l - x 3 ,c), {{1,3}, {2,4}}) 

+ B((xx, 1 - x 2 , 1 - x 3 ,c), (1 - x 1 ,x 2 ,x 3) 1 - c), {{1,4}, {2,3}}) 
< F'(xi 7 x 2 , l-x 3 ,l- c)F'(l -xi,l- x 2 ,x 3l c) 

+ F'(xi, 1 - x 2 ,x 3 , 1 - c)F'(l - Xi,x 2 , 1 - x 3 ,c) 

+ F'(xi, l-x 2 ,l- x 3 ,c)F'(l - Xi,x 2 ,x 3 , 1 - c) 
= F(xi,x 2 , 1 - x 3 )F(l - x-l, 1 - £2,2:3) 

+ F(a;i, 1 - x 2 , 2:3)^(1 - xi,x 2 , 1 - x 3 ) 

+ F(xu 1 — x 2 , 1 - x 3 )F(l - X!,x 2 ,x 3 ) 
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(2 =>■ 3) We can assume that F is not identically zero (otherwise, take two vertices 
of fugacity 0, and attach four outgoing edges to one of them - the isolated vertex can never 
be matched). Pick x e {0, l} 4 with F'(x\, X2, £3, X4) > 0. Lemma I3U1 implies that the flip 
F" of F' by x (1, 1, 1, 0) has a matchings circuit. By subdividing the i'th outgoing edge 
for each i with Xi = 1, we get a matchings circuit for F'. By Lemma [26] we get a matchings 
circuit whose weight-function is F. 

(3 => 1) Lemma [2^1 □ 

In particular by Theorem [S] and Proposition [5l Holant({R}) <ap #PerfMatch where 
R = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (1, 0, 1), (0, 1, 1)}. 
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